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We report new measurements of the neutron-deuteron elastic scattering cross section at energies 
from 135 to 250 MeV and center-of-mass angles from 80° to 130°. Cross sections for neutron-proton 
elastic scattering were also measured with the same experimental setup for normalization purposes. 
Our nd cross section results are compared with predictions based on Faddeev calculations including 
three-nucleon forces, and with cross sections measured with charged particle and neutron beams at 
comparable energies. 
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I. INTRODUCTION 

Three-nucleon systems comprise an area of long- 
standing interest in nuclear physics. The study of these 
systems has been enhanced by the growing database of 
precise measurements on two-nucleon systems [T], and 
the ability of modern potential models to provide ac- 
curate predictions of nuclcon-nucleon scattering observ- 
ables 2-5 . Furthermore, modern computational tech- 
niques have made it possible to calculate scattering cross 
sections and spin observables in three-nucleon systems 
for any kinematical configuration using the Faddeev for- 
malism [BJ, allowing the identification of experiments 
with strong sensitivity to the effects of three-nucleon 
forces (3NF). A comprehensive review of the theoreti- 
cal and experimental status with regard to 3NF in few- 
nucleon systems has recently been published by Kalantar- 
Nayestanaki, Epelbaum, Messchendorp, and Nogga [7]. 

In 1998, Witala and co-workers [5] employed Faddeev 
techniques to compute differential cross sections for nd 
elastic scattering at 65, 140, and 200 MeV, using solely 
two-nucleon forces as well as with the inclusion of the 
Tucson-Melbourne (TM) three-nucleon force [9. |. The 
cross section is seen to fall steeply with angle, due to the 
direct term, and then to rise again at backward angles 
due to the exchange term in the two-nucleon interaction. 
The 3NF term alone is nearly isotropic and thus is of 
magnitude comparable to or greater than that of the two- 
nucleon term for angles at which the cross section exhibits 
a minimum. Similar results were achieved by Deltuva, 



* |erol.ertan@boun.edu.tr| 

t Present Address: Aurora Flight Sciences, Cambridge, MA 02142 
i Present Address: Passport Systems Inc., Billerica, MA 01862 
§ Present Address: 9 Wagon Meadow Road, Santa Fe, NM 87505 



Machleidt and Sauer [5 in which a coupling of nucleon- 
A-isobar states to two-nucleon states is considered. Their 
study demonstrated that this coupling gives rise to an ef- 
fective 3NF force and a pion exchange mechanism. This 
suggests that measurements of the intermediate-to-large 
angle Nd cross section could reveal effects of 3NF. More- 
over, the 3NF predictions showed significant variation 
with incident energy indicating that energy-dependent 
measurements would be valuable. 

In recent years, several measurements have been car- 
ried out in this kinematic region, supplementing early 
measurements of nd cross sections near 150 MeV [10] 
and pd cross sections near both 140 MeV [11] [12] and 
200 MeV [HI H3] . 

Measurements of differential cross sections for dp elas- 
tic scattering at E d = 270 MeV performed at the RIKEN 
Accelerator Research Facility were reported in 2000 by 
Sakai et al. [H] and in 2002 by Sekiguchi et al. [15] . 
Measurements of pd elastic scattering cross sections have 
been performed at six energies between 108 and 190 MeV 
at the KVI facility [HJ |T7]. The cross sections at 
E p = 135 MeV were found to be 10-40% larger than 
those measured at RIKEN. Subsequent measurements 
by Sekiguchi el al. [18] at 135 MeV/ A using both pro- 
ton and deuteron beams support the original RIKEN 
measurements and contradict those from KVI. In 2008, 
Ramazani-Moghaddam-Arani et al. [TH] reported a new 
measurement of the pd cross section at KVI which yielded 
a result intermediate between the earlier KVI and the 
RIKEN cross sections. 

There have been fewer studies of the nd cross section 
than of the pd cross section. Mermod et al. (20] [21] mea- 
sured the differential cross section for nd scattering at 
95 MeV, and Maeda et al. [33] performed a measurement 
for center-of-mass angles from 10° to 180° at E n = 248 
MeV. These data were found to agree well with the pre- 
vious pd cross section measurements at E p = 252 MeV 
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reported in 2002 by Hatanaka et al. [2"3"j . 

Many of the previous experiments have been performed 
with polarized beams, and analyzing powers and other 
spin observables have been reported [HI [T7[ [TO] |2"2"H2"9"] . 
Discussion of these, along with measurements of the in- 
elastic Nd (breakup) cross section [30H33] . are beyond 
the scope of this paper. 

As stated earlier, most of the previous experiments 
have been carried out with charged particle beams, neces- 
sitating the consideration of Coulomb effects on the cross 
section. Also, with the exception of the KVI work [15111?] 
and the early work of Igo et al. [12] . previous measure- 
ments have been performed at a single energy. Moreover, 
there are still some lingering uncertainties in the magni- 
tude and shape of the differential cross section at 135 
MeV [19]. 

These factors have motivated the present measure- 
ment: a study of the ncutron-dcuteron clastic scatter- 
ing cross section, at large angles where the sensitivity to 
3NF is greatest, over a broad range of incident neutron 
energies. 

Since the absolute normalization of neutron scattering 
experiments can be a difficult problem, the np scattering 
cross section was also measured using the same experi- 
mental setup. The np data were then used to normalize 
the measured nd cross sections. Other recent nd mea- 
surements have used this technique to achieve precision 
sufficient to distinguish among calculations which display 
the explicit effects of 3NF [21]. 



II. EXPERIMENT 

The experiment was carried out at the Los Alamos 
Neutron Science Center (LANSCE) at the Los Alamos 
National Laboratory, Los Alamos, New Mexico. Neu- 
trons were produced as spallation products from an 800 
MeV H~ beam incident on a bare tungsten target. H~ 
pulses from the linear accelerator had a width of ap- 
proximately 0.2 ns, with an inter-pulse separation of 
1.8 us. This pulse structure allowed the energy of out- 
going (and incident) neutrons to be measured by time- 
of-flight techniques. The spallation method produced a 
"white" source of neutrons ranging from very low ener- 
gies up to nearly 800 MeV, with an energy distribution 
that depends on the production angle. 

This experiment was performed on a forward- angle, 
15° with respect to the primary H~ beam, flight path, 
to maximize the flux of neutrons in the energy regime 
above 100 MeV. The neutron beam was defined by two 
sets (horizontal and vertical) of Cu shutters with typical 
aperture of 3.8 x 3.8 cm, followed by a sweep magnet to 
filter out charged particles. The beam entered an evac- 
uated pipe containing a 2.7 m long steel collimator with 
an 1.3-cm diameter circular aperture surrounded by mag- 
netite shielding, and exited to pass through a 238 U foil 
fission ionization chamber |34| that monitored the beam 
flux as a function of neutron energy. Approximately one 



meter downstream from the fission chamber, the beam 
impinged on a cryogenic target cell containing either liq- 
uid deuterium (LD 2 ) or liquid hydrogen (LH 2 ). The 
target was located approximately 18 m from the spal- 
lation source. The intensity profile and position of the 
beam were measured at the target position by exposing 
a storage-phosphor image plate. 

The target geometry and orientation were designed 
to minimize the energy loss and reactions of outgoing 
charged particles. The target cell consisted of a hori- 
zontal cylindrical disk of 1.3 cm thickness and 12.7 cm 
diameter, with 51 /im Aramica entrance and exit win- 
dows. It was placed at an angle of 50° with respect to 
the incident neutron beam to provide an optimum path 
for the outgoing charged particles and to eliminate me- 
chanical constraints for neutrons exiting the cell at large 
angles on the opposite side of the beam. The target, 
located in an evacuated cylindrical scattering chamber 
with a 127 fj,m Kapton window, was cooled using a cryo- 
genic refrigeration system which employed gaseous He 
as its working fluid. The system had a nominal cooling 
capacity of 10 W at 20 K; a resistive heater was used to 
maintain a nominal absolute target pressure of 97 kPa. 

Scattered neutrons and recoiling charged particles were 
observed in coincidence. Protons and deuterons were de- 
tected by five telescopes, each consisting of a thin AE 
plastic scintillator backed by a pure Csl calorimeter. 
These detectors were positioned with their front faces 
100 cm from the center of the target, at mean laboratory 
angles of lab = 24°, 30°, 36°, 42° and 48°. The AE scin- 
tillators were 0.25 cm in thickness, and provided accurate 
information on the particle arrival time with an efficiency 
close to 100%. Their active areas of 10 x 10 cm 2 defined 
the solid angle for charged particle detection. The Csl 
calorimeters were 30 cm in depth, and provided a mea- 
sure of both the particle energy and its arrival time. 

The scattered neutrons were detected with five plas- 
tic scintillator bars, each 10 cm high x 10 cm thick x 
200 cm wide. The bars were stacked vertically to form a 
"wall" 200 cm wide x 50 cm high. Photomultiplier tubes 
were attached to the ends of each bar to allow both the 
neutron hit position and its time-of-flight to be deter- 
mined. The center of the neutron wall was positioned 
at a distance of 1.32 m from the target, and spanned a 
laboratory angle range from 34° to 108°. The face of the 
neutron wall was covered by four thin plastic veto scintil- 
lators to eliminate events produced by charged particles. 
The energy threshold of each bar was determined with 
low energy gamma-ray sources, and the photomultiplier 
gains were continuously monitored using cosmic-ray trig- 
gers. The experimental setup is illustrated in Figure [T] 
For nd elastic scattering, the average neutron angles in 
the laboratory and center of mass systems corresponding 
to each recoil deuteron detector are given in Table [T] 

Data acquisition was performed using standard NIM 
and CAMAC modules. A pre-scaled fraction of the "sin- 
gles" counts in each detector arm was read out by the 
electronics in addition to the coincidence events. Empty- 
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FIG. 1. Schematic diagram of the experimental setup. 



TABLE I. The average angles in lab and center-of-mass 
frames. 
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target runs were interspersed throughout the experiment 
to provide a measure of background. The target was 
filled with LD2 for the nd elastic scattering cross section 
measurements and with LH2 for normalization purposes. 
The experiment has been described in detail in |35j . 

III. DATA ANALYSIS 

Identification of np and nd elastic scattering events 
was achieved using a succession of cuts based on the 
charged particle and neutron detector pulse height and 
time information. The sequence began with a cut on the 
neutron beam energy. When the proton beam impinged 
on the tungsten spallation target, a time reference sig- 
nal, called t Q , was generated. Since gamma rays as well 
as neutrons and charged particles were produced in the 
spallation target, the "gamma flash" was easily identi- 
fiable as the leading peak in the time-of-fiight spectra 
such as that shown in Figure [2j The position of this 
peak could be used to establish the time offset calibra- 



tion for each detector with an uncertainty of ±0.67 ns 
(stat) ±0.20 ns (sys). The neutron beam energy for each 
event trigger was then deduced using the measured time 
between to and the trigger observed in the AE scintil- 
lator. This calculation assumes elastic scattering kine- 
matics, and includes a small time correction determined 
by Monte Carlo simulation for the charged particle en- 
ergy loss. The uncertainty in incident energy using this 
method was estimated to be 2 MeV. Scattering cross sec- 
tions were ultimately extracted in bins of 10 MeV width. 
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FIG. 2. Typical time-of-fiight spectrum for a AE detector, 
measured with respect to the arrival time of the proton burst 
at the spallation target (to). The "gamma flash" is observed 
at 60 ns. Charged particles are seen between 70 and 250 ns. 
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A. Neutron-proton Scattering Analysis 

Elastic neutron-proton scattering, using a LH2 target, 
was observed primarily in order to determine the neu- 
tron beam flux and the target thickness. Although both 
of these quantities were measured, several factors con- 
tributed to the uncertainty in these measurements. The 
fission chamber was well calibrated for neutron energies 
below 100 MeV, but uncertainties in the 238 U fission cross 
section limited the accuracy of the calibration at higher 
energies. The physical thickness of the cryogenic tar- 
get could be measured precisely at atmospheric pressure 
and room temperature, but not when the target was un- 
der vacuum and filled with liquid hydrogen or deuterium. 
Visual inspection of the target through the window of the 
scattering chamber under operating conditions revealed 
two effects, a bulging-outward of the cell windows and 
a steady stream of rising bubbles, both of which would 
change the effective target thickness. 

After the cut on the neutron beam energy, graphical 
cuts were applied on the AE — E histograms for the five 
charged particle telescopes. An example is shown in Fig- 
ure [3| Since the event rate was dominated by np elastic 
scattering, the recoil protons clearly stand out over back- 
grounds from protons elastically and inelastically scat- 
tered by other materials (curved band) and from protons 
experiencing strong interactions in the detectors as well 
as particles other than protons (horizontal band). For 
the data taken with hydrogen, these cuts were sufficiently 
selective to provide a good measure of the efficiency for 
charged particle detection, even when using the proton 
singles trigger in which information from the neutron de- 
tectors was not considered. This was particularly true 
after the backgrounds, determined from data taken with 
the target evacuated, were subtracted. 

The measured np scattering cross sections, with a 
"floating" absolute normalization, were compared to pre- 
dictions from the SAID multi-energy partial wave anal- 
ysis [I], known to be accurate at the 1% level in this 
energy regime based on comprehensive fits to the SAID 
NN database. Use of alternative partial wave analyses 
produced little difference in the results. Target thick- 
ness being the dominant factor in this renormalization, 
no strong energy dependence of the renormalization fac- 
tor is expected; hence an overall renormalization factor of 
1.25 was determined and was subsequently applied to the 
nd scattering data. Figure [4] shows the np cross section 
obtained from this work before and after renormalization. 

Assuming the fission chamber calibration to be sub- 
stantially correct, the correction factor represents a 
~25% increase in the effective target thicknessF] 



1 With regard to the validity of applying the same renormaliza- 
tion factor to the nd data: the fission chamber calibration and 
the window-bulging will certainly be the same for both hydro- 
gen and deuterium. Following the experiment, the target cell 
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FIG. 3. AE plotted versus E for np scattering at incident neu- 
tron energy 170 MeV and proton recoil angle 36°. The peak 
in the correlated (marked) band due to protons is clearly vis- 
ible. Protons of other energies can arise from elastic (higher 
energy) and inelastic (lower energy) scattering from, e.g., car- 
bon and oxygen nuclei in the windows of the target and scat- 
tering chamber. The horizontal band represents elastically 
scattered protons that have suffered strong-interaction energy 
losses in the Csl. These are included as valid events. 
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FIG. 4. The np cross section results with and without the 
correction factor of 1.25 ± 0.06. 



B. Neutron Detector Efficiency 

In addition to providing the normalization, the np 
scattering data were used to determine the efficiency of 
the neutron detectors. This is necessary since the nd 



was observed (at room temperature and atmospheric pressure) 
to have a permanent deformation; its physical thickness had in- 
creased by a factor of 1.45. The reduction of the effective target 
thickness due to bubbling apparently leads to the factor 1.25. 
In principle the reduction could be different for LH2 and LD2, 
since the boiling point temperature for LD2 is lower. However, 
the systematic error associated with the assumption of a constant 
renormalization factor is estimated to be 10% at worst, and has 
been included in the analysis. 
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data analysis requires the neutron-charged particle co- 
incidence trigger. Scattered neutron energies were de- 
termined by the measured beam energy and the proton 
scattering angle. The ratio of the number of coincident 
neutrons observed in the relevant regions of the wall of 
neutron bars to the total number of measured protons 
was used to determine the neutron detector efficiency. 
The results are presented in Figure [5] as a function of 
neutron energy. Also shown in the figure are the results 
of a Monte Carlo simulation [36]. The band shows the 
range of efficiencies which were found by fitting the data 
with a decaying exponential. Due to the effect of light 
attenuation in the scintillator, the efficiency decreases 
for neutrons incident near the ends of the bar. To ac- 
count for this position dependence, a separate analysis 
was performed in which the efficiencies for neutrons with 
the same energy but different positions were determined. 
These results provided correction factors for the data ob- 
tained at charged particle angles of 24° and 30°. The 
additional systematic error introduced by this procedure 
will be discussed later. 
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FIG. 5. Neutron detection efficiency in the scintillator bars as 
a function of neutron kinetic energy. The triangles represent 
the data obtained in the present measurement, the squares a 
Monte Carlo simulation [36] . The band is a fit to the data as 
described in the text. 



C. Charged Particle Detection Efficiency 

Recoiling charged particles will lose varying amounts 
of energy in the target depending on the location of the 
scattering event, and then will suffer further energy losses 
in the windows, air, and AE scintillator before being de- 
tected in the Csl. For nd scattering, a Monte Carlo simu- 
lation of these effects showed that all deuterons produced 
with energies greater than 35 MeV are detected. Since 
the lowest recoil deuteron energy in the present work is 
about 40 MeV, the efficiency of charged particle detection 
is assumed to be unity. 



D. Neutron-deuteron Scattering Analysis 

Analysis of the nd data was more challenging than that 
of the np data, since the nd elastic scattering cross section 
is much smaller than the d(n, np)n quasi-elastic scatter- 
ing cross section. Although in principle detection of a 
deuteron (in singles modc)in a charged-particle telescope 
is a unique signature of elastic scattering, in practice the 
limited resolution of these detectors does not allow un- 
ambiguous separation of deuterons and protons. This 
difficulty was resolved by employing the neutron-charged 
particle coincidence trigger. In coincidence mode, suc- 
cessive cuts were applied on both the difference between 
and the sum of the arrival times of the pulses at the left 
and right ends of the neutron bars. The difference pro- 
vides the location on the bar and the sum the total time- 
of-flight of the neutrons. The elastically scattered neu- 
trons associated with a given recoil deuteron angle satisfy 
kinematic conditions with respect to both position and 
time-of-flight, whereas the values for background events 
are distributed over the possible range. Typical AE — E 
histograms both before and after these conditions were 
imposed are shown in Figure [6] The suppression of the 
quasi-elastic events now allows the recoil deuteron peak 
to be observed. 

A typical measured distribution of neutron scattering 
angles for neutrons in coincidence with deuterons satis- 
fying the particle-identification cut is shown in Figure [7] 
In this histogram the neutron scattering angle is rep- 
resented by the measured difference between the pulse 
arrival times at the left and right ends of the scintilla- 
tor bars. Neutrons from the elastic nd events are eas- 
ily recognizable above the background. The region of 
interest is selected, and the events in that region are 
counted. For the majority of cases the background is 
relatively constant, allowing it to be easily interpolated 
under the region of the elastic peak and then subtracted 
from the total peak yield. In a few cases the background 
was not uniform, due to the residual quasi-elastic pro- 
tons which remained inside the deuteron cut. For these 
cases, a mathematical model was used to estimate the 
background. This procedure introduced additional sys- 
tematic errors, estimated to be around 10%. 

A correction factor was applied to each measurement 
to account for a slight degree of polarization in the neu- 
tron beam. The neutron polarization on the 15R flight 
path was measured previously with a CH2 polarimc- 
ter using left-right symmetric detection of np scattering 
yields [37] . The average vertical component of the beam 
polarization was found to be about 0.06, and to vary 
slightly with energy. The correction factor to the mea- 
sured yield is given by 1— p n *A y , where p n represents the 
polarization component normal to the scattering plane, 
and A y represents the analyzing power for the reaction. 
This correction was applied to both the np and the nd 
scattering data. The same partial wave analyses used 
in the evaluation of np cross sections were also used to 
determine np elastic analyzing powers to a precision of 
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FIG. 6. AE — E histograms of nd events for deuteron recoil 
angle 6d = 36° and incident neutron energy Tboam = 210 MeV 
a) before, b) after the cuts on the neutron timing as described 
in the text. 



the end, the corrections due to polarization are about 1%, 
and are not significant compared with the other errors. 



E. Systematic Errors 

An error of about 3% was assigned due to small devia- 
tions of the detectors from their assumed positions which 
was simulated with a Monte Carlo study. An error of 5% 
was due to np normalization. As mentioned earlier, the 
normalization coefficient may not be the same for LH 2 
and LD2 target due to different boiling points and bub- 
bling levels. Therefore, a generous error of 10% was in- 
troduced as there is no way of measuring this effect. The 
neutron detection efficiency error is about 5% on the av- 
erage, varying slightly with energy. A separate analysis 
showed that the error due to the unknown uranium cross 
section has to be included in addition to the np normal- 
ization error. The systematic error of interpolation in the 
unknown cross section region is estimated to be 9%. 

Additional errors were introduced for the most 
backward-angle data. The position dependence in the 
neutron detection efficiency led to systematic uncertain- 
ties of 14% at 24°, and 11% at 30°. Moreover, the method 
used to deal with nonuniform background in the neutron 
time-of flight difference spectrum at 24° introduces about 
10% systematic error. Therefore, the results at 36°, 42° 
and 48° carry about 16% total systematic uncertainty, 
whereas the results at 24° and 30° have about 24% and 
19% total systematic error, respectively. 



IV. RESULTS 
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FIG. 7. The neutron-bar time difference histogram at 9d 
36° and T bcam = 170 MeV. 



about 0.01. The estimated nd elastic analyzing powers 
had a greater uncertainty, since they were based on a 
limited set of available data. Because of the small size of 
p n , this uncertainty in A y constitutes a relatively small 
systematic uncertainty, even if no correction is made. In 



Cross sections for elastic nd scattering were obtained 
over the range of neutron (lab) energies from 130 to 255 
MeV. At lower energies, the energy loss of the recoil 
deuterons in the target and air was found to be too large; 
at energies above 255 MeV, particle identification based 
on the AE-E plot became unreliable. The data were 
binned into ± 5 MeV intervals at each of a total of eight 
centroid beam energies, and the differential cross sections 
were then calculated for each of the recoil deuteron an- 
gles. These results, and their uncertainties, are listed in 
Table [TTJ and are shown as angular distributions in Fig- 
ure [HJ The broad energy coverage apparent in these fig- 
ures is a unique feature of the present experiment, while 
the angular range of these new data is observed to fall in 
all cases within the region of the minimum of the cross 
section - i.e., where the sensitivity to the effects of three- 
nucleon forces is expected to be largest. 

Also included in these figures are the relevant data 
from previous nd, pd, and dp scattering experiments. The 
most extensive comparison that can be made with pre- 
vious work is with the KVI pd measurements of refer- 
ences [IHl [T7] . In nearly all cases these cross sections are 
larger than the present results. At our lowest energy, 135 
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MeV, the current results are in good agreement with the 
RIKEN dp cross sections [T4] . 

Figure [8] also displays the calculated cross sections. 
At 135 and 200 MeV, the Faddeev calculations of refer- 
ence [3] are shown both with and without the inclusion 
of 3NF. As expected, these results are significantly dif- 
ferent from each other only at angles larger than about 
90°. At both of these energies, the present cross sections 
are found to be in better agreement with the 3NF re- 
sults. At 150, 170, and 190 MeV the data are compared 
with calculations based on the CD Bonn potential with- 
out and with the inclusion of a A-N component [5] . At 
the largest angle where this effect is most significant, the 
data favor its inclusion. 



V. CONCLUSION 

In this study, the differential cross section for nd elastic 
scattering was measured in a continuous incident neutron 
energy range from 135 to 250 MeV, by detecting scat- 
tered neutrons and recoil deuterons in coincidence, with 
the aim of elucidating the contribution of three-nucleon 
forces (3NF), in particular the energy dependence of this 
effect. The absolute scale of the nd cross section was de- 



termined by concurrent measurements of np elastic scat- 
tering, which were also used to determine the neutron 
detection efficiency. 

The effect of 3NF is clearly seen in this work. The data 
at angles near the minimum in the cross section, where 
the 3NF contribution is most effective, are in excellent 
agreement at all energies with the theoretical predictions. 

The 3NF effect could be further tested by confronting 
the present data with theoretical predictions for the dif- 
ferential cross section at fixed angles as a function of 
incident neutron energy. 
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